The HeH^ molecular ion in a magnetic field. 

A. V. Turbiner* and N. L. Guevara^ 
Instituto de Ciencias Nucleares, Universidad Nacional Autdnoma de Mexico, 
Apartado Postal 70-543, 04510 Mexico, D.F., Mexico 
(Dated: February 2, 2008) 

Abstract 

A detailed study of the low-lying electronic states ^S,^S,^n, of the HeH"^ molecular ion in 
parallel to a magnetic field configuration (when a-particle and proton are situated on the same 
magnetic line) is carried out for B = — 4.414 x 10^'^ G in the Born-Oppenheimer approximation. 
The variational method is employed using a physically adequate trial function. It is shown that 
the parallel configuration is stable with respect to small deviations for S-states. The quantum 
numbers of the ground state depend on the magnetic field strength. The ground state evolves from 
the spin-singlet state for small magnetic fields B < 0.5 a.u. to the spin-triplet unbound state 
for intermediate fields and to the spin-triplet strongly bound ^11 state for B > 15 a.u. When the 
HeH"*" molecular ion exists, it is stable with respect to a dissociation. 
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I. INTRODUCTION 



It is well-known that a very strong magnetic field may appear on the surface of neutron 
stars, B ^ 10^^ — 10^^ G. It seems natural to expect that in a vicinity of the surface some 
atomic and molecular systems can occur. For many years this assumption has motivated a 
development of atomic-molecular physics in a strong magnetic field hoping that it might lead 
to understanding of the spectra of these compact stars and be important for a construction 
of a model of the atmosphere. Recent revolutionary observation of the absorption features 
in a spectra of the isolated neutron star 1E1207. 4-5209 by Chandra X-ray observatory [1] 
brought an extra attention to the field. 

So far, a major attention has been paid to a study of one-electron atomic and molecular 
systems (for a review, see [2]). During the last years the existence of several exotic, strongly 
bound one-electron molecular systems has been predicted. In particular, it was shown that 
for a sufficiently strong magnetic field B > 5 x lO-'^^G the exotic ion (HeH)^~^ can exist 
in parallel configuration (the magnetic field is directed along internuclear axis) as optimal 
as well as its excited states Itt, IS [3]. A complete list of possible one-electron atomic and 
molecular systems is known for B < Bgchwinger — 4.414 x 10^^ G. In turn, the exploration of 
two-electron systems has been mostly restricted to a study of atomic- type systems H~, He 
with the only exception of the H2 molecule (see [4], [5] and references therein). Recently, 
a first detailed study of the molecular ions H3 [6] and [7] has been carried out (for a 

brief review of a general situation sec [8]). It was shown that for large magnetic fields if such 
an ion exists its optimal configuration is the linear parallel one - heavy charged particles are 
situated along a magnetic line. A goal of the present paper is to make a study of the Coulomb 
system (apee) in a magnetic field ranging in B — — 4.414 x 10^^ G and thus to estabhsh the 
possible existence of the simplest two-electron diatomic heteronuclear molecular ion HeH''" 
and to explore its lowest excited states. In the field-free case the theoretical existence of 
HeH^ was established long ago by Beach [9] (sec also [10]) while its several excited states 
were found (see, e.g., [11]). Experimentally, this molecular ion was first observed in the mass 
spectra of helium-hydrogen mixtures (see, e.g., [12]). The HeH"*" ion is thought to be among 
the first molecular species to appear in the Universe (see [13] and references therein). An 
attempt to explore the HeH''" ion in ultra-strong magnetic field was made in [14] . 

Atomic units are used throughout {h—rrie—e—l), although energies are expressed in Ry- 
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dbergs (Ry). The magnetic field B is given in a.u. with Bq — 2.35 x 10^ G. 



II. GENERALITIES 

We study the system (apee) placed in a uniform constant magnetic field. It is assumed 
that the a-particlc and proton are infinitely massive (Born-Oppenheimer approximation of 
zero order). It seems natural from physical point of view that the optimal configuration is 
parallel when the a-particle and proton are situated along the magnetic field line. Never- 
theless a stability of this configuration with respect to small perturbations will be checked. 

The Hamiltonian which describes the system (apee) when the magnetic field is oriented 
along the z direction, B = (0, 0, B) is 

^-E(P^ + -^^)'- E ^ + - + | + 2B.S, (1) 

K=a,p 

(see Fig. 1 for the geometrical setting and notations), where = —iVe is the momentum 
of the i?-th-electron, the index k runs over the a-particle and proton, Za = 2, Zp = 1 and 
ri2 is the interelectronic distance, S = Si + §2 is the operator of the total spin. The vector 
potential Ae which corresponds to the constant uniform magnetic field B is chosen in the 
symmetric gauge, 

1 B 

Ae = -{Bx r,) = -(-?/,, xe, 0) . (2) 
Finally, the Hamiltonian can be written as 

i=l ^ ^ e=i,2 ^'^ '■^ 

K=a,p 

where — Izi + lz2 Sz — Sz^ + Sz2 are the z-components of the total angular momentum 
and total spin, respectively, and p| = x| + , 1—1,2. 

The problem is characterized by two conserved quantities: (i) the operator of the z- 
component of the total angular momentum (projection of the angular momentum on the 
magnetic field direction) giving rise to the magnetic quantum number m and (ii) the operator 
of the 2;-component of the total spin (projection of the total spin on the magnetic field 
direction) giving rise to the spin quantum number m^. Hence, any eigenstate has two 
quantum numbers assigned: the magnetic quantum number m and the spin quantum number 
rris. 
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FIG. 1: Geometrical setting for the HeH"*" ion parallel to a magnetic field directed along the z-axis. 
The a-particle and the proton (marked by bullets) are situated on the z-line at distance R. 



As the magnetic field increases a contribution from the linear Zeeman term (interaction 
of the spin with a magnetic field B • S) becomes more and more important. Hence, it is 

natural to assume that for large magnetic fields the states with minimal (when both 
electron spins are antiparallel to the magnetic field) will represent the states with lowest 
total energies. 

The space of eigenstates is split into subspaces (sectors) with each of them characterized 
by definite values of m and m^. Thus, to classify eigenstates we follow the convention 
widely accepted in molecular physics using the quantum numbers m and m^. The notation 
is ^^'^^M, where 25' + 1 is the spin multiplicity and corresponds to 1 for singlet 5 = and 3 
for triplet 5 = 1, as for the labeling M it is used the Greek letters S, 11, A that correspond 
to the states with \m\ = 0, 1, 2, respectively. We study only states with m — 0, —1, —2. 

As a method to explore the problem we use the variational procedure. The recipe of 
choice of trial functions is based on physical arguments [15] . A trial function for states with 
orbital magnetic quantum number m is chosen in the form 

^{trial) _ ^-^ ^^p^^-^p\M^irn(t>i^-airiA-a2riB-a3r2A-0!4r2B-B/3ipy4:-B/32P^/4:+'yri2 ^^-^ 

(cf.[4]), where (Te = 1, —1 stand for spin singlet {S — 0) and triplet states {S — 1), respec- 
tively. The Pi2 is the operator which interchanges electrons (1 <-> 2). The parameters q;i_4. 
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(3i-2 and 7 as well as the internuclear distance R are eight variational parameters. 

Calculations were performed using the minimization package MINUIT from CERN-LIB. 
Multidimensional integration was carried out using a dynamical partitioning procedure: a 
domain of integration was divided into subdomains following the integrand profile and then 
each subdomain was integrated separately (for details see, e.g., [2]). Numerical integration 
was done with a relative accuracy of ~ 10~^ — 10~^ by use of the adaptive DOIFCF routine 
from NAG-LIB. Computations were performed on a dual DELL PC with two Xeon processors 
of 2.8 GHz each (ICN) and dual DELL PC with two Xeon processors of 3.06 GHz each 
(CINVESTAV). We implemented a parallelization in the routine for the calculation of the 
integrals. Minimization process for each particular magnetic field and for a given state was 
quite painful and lengthy, it took from 100 to 1000 hours of CPU time, however, when a 
minimum is found it was needed a few minutes to compute a minimal variational energy. 

III. RESULTS 

In absence of a magnetic field the HeH"^ molecular ion exists and it is stable. Its ground 
state is the ^T, {S = and m = 0) state [9]. In order to establish the existence of HeH"*" in 
a magnetic field wc carried out a detailed study of the Coulomb system (apee) in parallel 
configuration in the state for magnetic fields B = — 10000 a.u. (see Table I). For this 
case the variational trial function ■0*''*"' (4) is used at CTg = 1 and m — 0. 

The calculations indicate the optimal geometry is parallel (a-particle and proton stay on 
the same magnetic line). For this configuration for all studied magnetic fields a minimum 
in the total energy curve Et{R) is developed. Furthermore, we checked that this minimum 
is stable with respect to small inclinations of the internuclear axis from the magnetic field 
direction. It is found that with an increase of the magnetic field strength the total energy 
of the {apee) system in the state grows, it becomes more bound (in particular, the 
double ionization energy increases Ej — 2B — Et) and the system becomes more compact 
(the internuclear equihbrium distance decreases), see Table I. For B < 100 a.u. the {apee) 
system in the ^Ti state is stable with respect to the dissociation He+ + H or He + p. The 
vibrational and rotational energies are calculated for B = 0.5 a.u. The vibrational energy is 
larger than the rotational energy which is typical for field-free case. 

We carried out a detailed study of the spin-triplet state {S — 1 and m = 0) of the 
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TABLE I: HeH+ ion in the state S: Total Et, double ionization E[, vibrational Eq and rotational 
Eq"* energies in Ry, and equilibrium distance in a.u. of the (apee) system in the state ^E. Magnetic 
fields for which the state is the ground state are marked by * (see text). " Result from [9] 
(rounded). Total energy E{*) = £'(He"*'(ls)) + E(H(ls)) in Ry with electrons in spin-singlet state 
is from [2]. Total energy of He are from [16], ^ marks a result of linear interpolation. 



B (a.u) 


Et 


Ei 


Req 


Tpvib 


rprot 




£;(i^e(1^0+)) 


0* 


-5.943 


5.943 


1.451 


0.016 




-5.0 


-5.807 




-5.957 " 




1.463 










0.5* 


-5.843 


6.843 


1.430 


0.018 


0.005 


-4.8>^ 


-5.6>^ 


1 


-5.626 


7.626 


1.404 






-4.544 


-5.459 


10 


5.337 


14.663 


0.789 






7.725 


6.129 


100 


166.513 


33.487 


0.370 






173.327 


169.837 


1000 


1927.35 


72.65 


0.212 






1944.37 




10000 


19848.67 


151.30 


0.0928 






19893.62 





(apee) system in parallel configuration in the domain of magnetic fields < B < 10000 a.u. 
For this state the variational trial function (4) with ae — —1 and m = is used. 

Similar to ^T, all total energy curves Et{R) of the {apee) system for magnetic fields 
B — {)— 10000 a.u. in the state displays a minimum. However, to the contrary to the 
state the minimum is always shallow, see Fig. 2-5 as an illustration. It is found that with 
an increase of the magnetic field strength the total energy decreases. The system gets more 
"compact" - the internuclear equilibrium distance of a shallow minimum decreases slowly, 
see Table II. We expect that the system in the ^E state "exists" in a form of two separated 
atomic-hke systems He''"(ls) -|- H(ls) situated at large distances with electron spins directed 
against a magnetic field direction. 



6 



A comparison of the total energies of the (apee) system in and states in equihbrium 
indicates to the appearance of a level crossing (see Tables I and II). This crossing of the 
^E and ^E states occurs at B ki 0.5 a.u. It implies that the ■'^E state is the lowest energy 
state aX B < 0.5 a.u., where the (apee) system is stable towards possible decays. Hence, 
one can conclude that the stable HeH"^ molecular ion exists at S < 0.5 a.u with state 
as the ground state. For B > 0.5 a.u. the situation is changed - the ^E state may become 
the ground state of the HeH'*' molecular ion. In fact, the state ^E is at most metastable 
with very short lifetime or very likely is simply unbound - the HeH"*" ion dissociates to 
He'*"(ls) + H(ls) with the electrons in spin-triplet state. A general conclusion can be drawn 
that for 0.5 < B < 15 a.u. the system (apee) likely exists in a form of two separate atomic- 
type systems He+(ls) + H(ls) with total electron spin equal to one. 

As next we carried out a detailed study of the state ^11 (5 = 1 and m = —1) of the (apee) 
system in parallel configuration for magnetic fields < B < 4.414 x 10^^ G (see Figs. 2-5). 
For this state the same variational trial function ■0*^*"' (4) is used at (7e = — 1 and m — —1. 

The results of the calculations indicate unambiguously the existence of a deep minimum 
in the total energy curve Et{R) of the (apee) system in parallel configuration for all studied 
magnetic fields < B < 4.414 x 10^'^ G. Hence, it can be concluded that the HeH^ molecular 
ion may exist in the ^H state. Table III contains the results for the total energy Et and 
the internuclear equilibrium distance R^q for the state ^H. Perhaps, it is worth mentioning 
that the previously obtained results for ultra-magnetic fields [14] are less accurate for both 
studied two-electron systems: (i) the HeH"^ ion in the present calculation and (ii) the He 
atom in the calculation made in [17]. With an increase of the magnetic field strength the 
total energy decreases, the system becomes more bound (double ionization energy increases 
being equal for triplet states to Ej = —Et) and more compact (the internuclear equilibrium 
distance decreases). Table HI shows the total energies of three dissociation channels of 
HeH+ in =^H state: ^ He(l=^(-1)+) +p, ^ He+(ls) + H(2p_i) and ^ He+(2p_i) + H(ls) for 
different magnetic fields. It can be seen that the state ^H is stable towards any dissociation. 
It is worth noting that the dissociation energy needed increases monotonously in all range of 
studied magnetic fields. Therefore the stability of the system with a magnetic field increases. 

One can see that the crossing of the total energies of the '^H and states at equilibrium 
occurs at B 15 a.u. For larger magnetic fields the total energy of the ^H state is lower 
than the total energy of ^E. Hence, the ground state for the HeH"*" ion for B > 15 a.u. is 
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TABLE II: HeH^ ion in the state '^S: Total Et, double ionization Ei energies in Ry and equilibrium 
distance in a.u. of the {apee) system in the state ^S. Magnetic fields for which the state might 
be the ground state are marked by *. Total energy E{*) = £^(He"^(ls)) + E(H(ls)) (in Ry) for the 
electrons in spin-triplet state is from [2] . Total energies for the He atom with total electron spin 1 
are from [16]. " Rcf [11]. 



B(a.u) 


Ef 




E{*) 


E{He){l^O+) 





-5.0064 


4.50 


-5.0 


-4.3504 




-5.00878" 


4.35 






0.5* 


-5.8712 


4.32 


-5.8637 


-5.3796** 


1* 


-6.5524 


4.04 


-6.5443 


-5.3013 


10* 


-12.209 


2.70 


-12.275 


-9.255 


15* 


-13.936 


2.47 


-14.043 


-10.622** 


100 


-26.172 


1.68 


-26.673 


-19.686 


1000 


-54.020 


1.20 


-55.628 




10000 


-103.07 


0.83 


-106.38 





realized by the ^11 state. Eventually, one can draw a conclusion that the state is the 
lowest energy state of the {apee) system for 0.5 < B < 15 a.u. and the HeH"'' ion can exist. 

In order to complete the study we consider the state (5" = 1, m = — 2 ) of the {apee) 
system in the domain of magnetic fields < -B < 4.414 x 10^^ G. It is worth mentioning 
that the state '^A for the ficld-frcc case was studied in [11] and [18]. For this state the 
same variational trial function -0*^*"^ (4) is used but with Ue — —1 and m — —2. For 
all total energy curves Et{R) for different magnetic fields the minimum is found. Hence, 
the ion HeH''" may exist in ^A state for all magnetic fields. Table IV shows the total Ex 
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TABLE III: HeH^ ion in the state '^11: Total Et and vibrational Eyn, energies in Ry and equilibrium 
distance Req in a.u. of the HeH"*" ion in the state. Magnetic fields for which the state is 
the ground state are marked by *. Total energies E{**) = E(He^{2p-i) + H(ls)) and E{* * *) = 
£'(He"^(ls) + H(2p_i)) (in Ry) for the electrons in spin-triplet state from [2]. Total energy for 
He atom in spin-triplet electronic state in magnetic field is from [17] (rounded). * Result of 
interpolation. " from [11]. ^ from [16]. from [14]. 



B(a.u) 


Ej' 


Req 


Evil) 


E{**) 


E(* * *) 


E{He{lH-ir)) 





-4.263 


7.61 




-2.0 


-4.25 


-4.266^ 




-4.276'' 


7.7 










1 

J. 


^ Q7^ 










Q1Q 


10 








O. ( 


1 -L.UOU 




15* 


-13.96 


0.925 


0.04 




-12.66 


-13.58* 


100* 


-28.36 


0.440 


0.15 


-20.30 


-24.37 


-26.15 


1000* 


-64.24 


0.203 


0.50 


-43.98 


-51.60 


-56.06 




-63.07^= 










-54.95*= 


10000* 


-133.49 


0.104 


1.41 


-86.79 


-99.98 


-110.30 




-129.7^= 










-106.9^ 


4.414 X 10^^ G * 


-160.50 


0.092 











energy and the internuclear equilibrium distance Rgq for the state ^A. With an increase 
of the magnetic field strength the total energies decrease, the system becomes more bound 
(double ionization energies increase, Ej = —Et, for triplet states) and more compact (the 
internuclear equilibrium distance decreases). At 5 ~ 300 a.u. a crossing of the total energies 
of and at equihbrium occurs. For larger magnetic fields the total energy of the ^A is 
smaller than the state. It imphes that for B > 300 a.u. the ^ A state is the lowest excited 
state. 
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TABLE IV: Total energy in Ry and equilibrium distance in a.u. of the HeH"*" ion in the state ^A. 



B(a.u.) 






n 


-4 108 




1 

± 


O .\J\Jrt 


4 71 7 


10 


-11.02 


1.581 


100 


-25.78 


0.486 


1000 


-58.68 


0.218 


10000 


-123.62 


0.112 


4.414 X 10^^ G 


-149.0 


0.096 




IV. CONCLUSION 

The low-lying electronic states ^E,^E,^n,^A of the Coulomb system apee in parallel 
configuration are studied for 5 = — 4.414 x 10^^ G using the variational method in the 
Born-Oppcnheimer approximation of the zero order. It is shown that for all studied states 
and all magnetic fields the parallel configuration is always optimal being stable with respect 
to small inclinations. The state of the lowest total energy depends on the magnetic field 
strength. It evolves from the spin-singlet state for small magnetic fields B < 0.5 a.u. to 
spin-triplet ^E (unbound state) for B > 0.5 a.u. and, finally, for B > 15 a.u. to spin-triplet 
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1 

1 
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\ 






B=100 a.u. 


E^(Ry) 
-24 




n 








-28 

















0.5 




1 


R(a.u.) ■'■5 



FIG. 3: Total energy curves for triplet states and ^11 for a magnetic field B = 100 a.u. 
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-60 
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0.4 


0.8 


1.2 



R(a.u.) 

FIG. 4: Total energy curves for triplet states and ^11 for B = 1000 a.u. 

state (strongly bound state). This result shows that the HeH"*" molecular ion exists as a 
bound state for magnetic fields B < 0.5 a.u. and B >15 a.u. in a parallel configuration and 
it is stable. In a domain 15 < B < 300 a.u. the state having the total energy next after the 
one of the ground state is the unbound (repulsive) state. For B > 300 a.u. the ^ A state 
is the lowest excited state of the HeH"*" molecular ion. As for 0.5 < 5 < 15 a.u. the system 
{apee) likely exists in a form of two separate atomic-type systems He^(ls) + H(ls) with 
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-40 

E.^(Ry) 

-80 



-120 

0.2 0.4 0.6 0.8 

R(a.U.) 

FIG. 5: Total energy curves for triplet states and ^11 for B = 10000 a.u. 

total electron spin equal to one which are situated at very large distance from each other. 
It is worth emphasizing that for B > 2000 a.u. both one- and two-electron ions {HeHY'^ 
and (HeH)^ can exist. As for the ion (HeH)'^^ it becomes stable at B > 10000 a.u. 
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